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Models of the platelet-derived growth factor (PDGF)-nke 
domains of vascular endothelial growth factor (V£GF) and 
placenta growth factor (PIGF) were built based on their 
homology to POGF. These domains contain most of the 
determinants for receptor binding. The sequences of these 
proteins exhibit limited but significant homology to that of 
platelet-derived grovyth factor (PDGF), a member of the 
cystine knot growth factor family. The eight cysteine 
residues that are involved In intra- and interchain disul- 
pbide bonds are conserved. Two high affinity receptors for 
VEGF have been identified, only one of which binds PIGF. 
The models presented here are consistent with results that 
show that VEGF receptor binding is mediated by charged 
residues in the loops. A comparison of the models suggests 
that the difference in receptor-binding specificity between 
VEGF and PIGF may be due to differences in the distribu- 
tion of positively charged residues and the exposure of 
hydrophobic residues in the loops. 
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Introduction 

In the early 1970s, it was hypothesized that tumour growth 
depends on angiogenesis and that the inhibition of angio- 
genesis might be an effective approach to the treatment of 
solid tumours (Folkman, 1971). There is now much evidence 
to support this proposal. Tumours become vascularized by 
inducing angiogenesis in nearby blood vessels, causing them 
to produce capillaries that grow into and vascularize the 
tumour. At least 20 peptides have been shown to be directly 
or indirectly angiogenic but there is mounting evidence that 
the primary regulator of tumour angiogenesis is vascular 
endothelial growth factor (VEGF) (Ferrara et ai, 1993). VEGF 
was originally described in highly vascularized tumours where 
its expression is increased by hypoxia (Shweiki et a/., 1992) 
and was first isolated from bovine pituitary folliculosCellate 
cells as a secreted heparin-binding protein of 46 kDa molecular 
weight (Ferrara and Henzel, 1989). It is a specific mitogen 
for endothelial cells and is a potent inducer of vascular 
permeability, which allows tumour cells to enter blood vessels 
(Ferrara and Henzel, 1989; Houck et oL, 1991). The growth 
of tumours in vivo can be suppressed by monoclonal antibodies 
raised against VEGF but these antibodies have no effect on 
the growth of tumour cells in vilw (Kim et ai, 1993). VEGF 
is overexpressed in glioblastoma cells and endothelial cells in 
tumours exhibit a 2-3-fold increased expression of the VEGF 
receptor in vivo (i'late et af.^ 1992). Glioblastoma growth 
in vivo can be inhibited by trans fecting endothelial cells with 



a dominant negative mutant of the VEGF receptor (Millauer 
et al., 1994). 

The molecular cloning of human cDNA for VEGF yielded 
four transcripts encoding species with insertions or deletions 
in the cDNA occurring at a common site (Tischer et aL, 1991). 
The four transcripts encoded a 26 amino acid hydrophobic 
signal sequence and had identical mature amino termini. The 
transcript encoding the 46 kDa form generates a 165 amino 
acid peptide following signal cleavage; this peptide is known 
as VEGF165. The other transcripts encode VEGF121, which has 
a 44 amino acid deletion between positions 116 and 159 of 
VEGF|65, VEGF|g9, which has 24 amino acids inserted at 
position 116, and VEGF206, which has 41 amino acids inserted 
at position 116 (Tischer et al, 1991). 

Two tyrosine kinases, fms-like tyrosine kinase (flt-I) and 
the kinase insert domain-containing receptor (KDR) have been 
identified as the high affinity receptors for VEGF (de Vries 
et al., 1992; Terman et al, 1992). The murine homologue of 
KDR is flk-1, whose sequence is 85% identical with diat of 
KDR. The fit-l and KDR genes encode tyrosine kinase 
receptors consisting of an extracellular domain containing 
seven immunoglobulin domains, a transmembrane domain 
and a split tyrosine kinase intracellular domain (Shibuya 
et al. 1990). 

The longer forms of VEGF bind heparin through a site 
that appears to be distinct from the receptor-binding domain 
(Vaisman et al., 1990). The heparin binding site is located in 
the basic C-terminal region of the longer forms of VEGF since 
VEGF|55 can be cleaved by plasmin to yield VEGF|m, which 
lacks residues 111-165 and has no affinity for heparin (Park 
et fl/., 1993). Since the sequence of residues l-llO displays 
limited but significant homology (19%) to that of PDGF, this 
suggests that the first 1 10 residues of VEGF fold independently 
to form a PDGF-like domain that binds to the VEGFR while 
the 55 C-terminal residues form the heparin-binding domain. 
VEGF 12 1 lacks this basic region and does not bind heparin but 
can trigger angiogenesis and vascular permeability (Houck 
et al., 1992; Kondo et al., 1995). 

The carboxy-terminal peptide (residues 111-165) produced 
by plasmin cleavage of VEGF155 has no affinity for soluble 
KDR or flt'l receptors (Kcyt et al„ 1996a). VEGF„o binds to 
soluble KDR with similar affinity to VEGF165 and VEGF,2i 
in the absence of heparin (Keyt et al., i996a). By contrast, 
VEGF,2j and VEGF|,o bind to fit- J witli 10-20-fold decreased 
affinity compared with VEGF165 (K.eyt et al, 1996a). A 
VEGFj65/VEGF||o heterodimer has an affinity for flt-l inter- 
mediate between VEGFuq and VEGFj^s homodimers (Keyt 
et al, 1996a). These results show that the C-terminal domain 
of VEGF contains determinants for binding to flt-I even though 
the 55 amino acid peptide by itself has no detectable affinity 
for ;?r-y. 

In addition to VEGF, placenta growth factor (PIGF) has 
been shown to bind to flt-1 with high affinity, but not to KDR 
(Kendall et al., 1994). PIGF induces procoagulant tissue factor 
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production by monocytes and endothelial cells and is a very 
weak stimulator of endothelial cell chcmotaxis and proliferation 
(Clauss et ai, 1996). The sequence of PIGF is 47% identical 
with that of VEGF. Alternate splicing of mRNA for PIGF 
generates two variants containing 132 and 153 amino acids 
(Maglione et aL, 1993). known as PIGF-1 and PIGF-2. 
respectively. The insertion of a 21 amino acid basic region 
near the C-terminus of the longer form confers heparin-binding 
ability on it (Hauser and Weich, 1993). 

The active forms of VEGF and PIGF are dimeric; in 
addition to VEGF and PIGF homodimers, a heterodimer of 
VEGF and PIGF has been purified and characterized (DiSalvo 
et aL, 1995). Heterodimcrs composed of combinations of 
VEGF,2i, VEGF, 55. PlGF-l and PiGF-2 have been synthe- 
sized in vitfv from monomers and shown to have mitogenlc 
activity (Birkenhager et al, 1996). That the various hetero- 
dimcrs are biologically active suggests that the determinants 
of dimerization are located in the first -110 residues of VEGF 
and PIGF. Homodimeric PIGF has very weak mitogenlc 
activity compared with homodimeric VEGF and the VEGF/ 
PIGF heterodimer (DiSalvo et ai, 1995); this is possibly 
because the KDRJflk-i receptor can transduce a mitogenic 
signal while activation of the /it- 1 receptor does not seem to 
produce a similar response (de Vries et al., 1992; Millauer 
et uL, 1994). 

The VEGF and PIGF genes display a limited sequence 
similarity to thai of platelet-derived growth factor (PDGF). a 
member of the cystine knot growth factor family that also 
includes transforming growth factor p (TGPP), nerve growth 
factor (NGF) and chorionic gonadotropin (CG) (McDonald 
and Hendrickson, 1993). The characteristic fold of this family 
is a pair of twisted anliparallel |3 sheets linked by a cystine 
knot formed by three disulphide bonds (Murray-Rust et al, 
1993). One of the disulphide bonds penetrates a ring formed 
by the other two bonds and the backbone atoms linlcing them. 
The sequence similarity in the family is low but the residues 
involved in the cystine knot are conserved. The structures are 
most similar about the knot and diverge rapidly from each 
other as the distance from the knot increases. The P strands 
have similar hydrogen-bonding patterns in each structure. The 
cystine knot growth factors are dimeric in vivo but the mode 
of dimerization is different in each case. The mechanism of 
action is thought to involve ligand-induced dimerization of the 
tyrosine kinase receptors. 

The cysteines involved in forming the cystine knot are 
conserved in the sequences of VEGF and PIGF. Two other 
cysteines diat are essential for dimerization of PDGF are also 
conserved in VEGF and PIGF. These cysteines have been 
shown to be essential for the dimerization and biological 
activity of VEGF (Polgens et aL, 1994). 

Because of the low sequence similarity among the cystine 
knot growth factors, a comparison of the sequences of VEGF 
and PIGF with those of PDGF-A and PDGF-B gives little 
insight into the receptor-binding determinants. Alanine- 
scanning mutagenesis has been used to identify charged 
residues involved in VEGF binding to flf-l and ICDR (Keyt 
et aL, 1996b). VEGF binding to K.DR/flk-1 is mediated by a 
positively charged surface while binding to fit- J is due to 
negatively charged residues. A model of VEGF based on the 
structure of PDGF shows that these residues are clustered in 
loops at opposite ends of the VEGF dimer (Keyt et aL, 1996b). 
This model was constructed by aligning the sequence of VEGF 
with that of PDGF and substituting the side chains of the 



charged residues involved in receptor binding for the corres- 
ponding residues in PDGF. 

In the present work, models of the PDGF-like domains of 
VEGF and PIGF have been built by homology modelling 
using the crystal lographic structure of PDGF-B B as a template. 
Deletions in the loop regions were modelled using database 
searching and a systematic classification of P hairpin structures 
(Sibanda et ai, 1989). The models were refined using energy 
minimization and molecular dynamics. The correctness of 
the models was assessed by evaluating the compatibility of 
the sequences with the cystine knot fold. This was done by 
comparing the models with their sequences using a 3D profile 
computed fi-om the atomic coordinates of the models and 
calculating the compatibility of each residue in the sequence 
with the environment defined by its position in the profile. 
The sum of the compatibility scores over all residues in the 
sequence gives the self-compatibility score; it has been shown 
that the 3D profiles of correct protein structures match their 
sequences with high self-compatibility scores (Luthy et aL, 
1992). 

Given the receptor-binding affinity of VEGF|,o, this model 
should contain all the determinants of K.DR binding and the 
most important determinants of flt-l binding. The model of 
VEGF is consistent with the available data on the residues 
that mediate receptor binding. No corresponding data for PIGF 
are available at the time of writing. It is hoped that the models 
will stimulate and guide further work to identify the residues 
that detennine receptor binding in VEGF and PIGF. The 
models suggest an explanation for the difference in receptor- 
binding specificity of VEGF and PIGF. 

Materials and methods 

The X-ray structure of PDGF-BB at 3,0 A resolution was used 
as the basis of the models (Oefner et aL, 1992). Monomer A 
of the PDGF structure was used as the template since this 
contains the more complete set of coordinates. Sequences were 
aligned using Clustal (Higgins and Sharp, 1988) and manually 
adjusted to preserve the P sheet structure (Figure 1 3). Molecular 
graphics and model building were done using the InsightH 
software package (Biosym Technologies, 1994) running on a 
Silicon Graphics Indy workstation. Attempts to predict the 
secondary structure of VEGF and PIGF using the PredictProtein 
(Rost, 1996) and SOPM algorithms (Geourjon and Deleage, 
1995) did not produce useful information. 

The PDGF(A) structure corresponds to residues 19-109 of 
VEGF/PIGF. There are no coordinates for the residues in 
PDGF(A) corresponding to residues 38^ of VEGF/PIGF. 
This is due to proteolytic cleavage in this region of PDGF, 
which makes it impossible to resolve the structure of loop 1. 
The Homology module of the Insightll package was used to 
build coordinates for residues 17-37, 45-83 and 89-109, using 
PDGF(A) as the basis. Loop 11 (residues 62^4) is the same 
length in PDGF and VEGF/PIGF, so this loop was built using 
the loop from PDGF. Loops 1 and III were built by database 
searching using the Homology module. Loop III (residues 84- 
88) is a P hairpin in the growth factors and was modelled 
using the published systematic classification of P hairpins 
(Sibanda et aL, 1989) as a guide. Each loop was chosen to 
have a minimum root mean square (r.m.s.) distance between 
the tails of the loop and the sheets and to have the maximum 
sequence similarity with VEGF/PIGF. In addition, the com- 
patibility of residues in the loop with the environments provided 
by the fragments was assessed using die Profiles program. The 
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Fig. I. MOLSCRIPT (Kraulis. 1991) stereo drawing of the VEGF dimer Hydrophobic residues in the dimer interface are drawn in ball-ond-stick 
representation. 



models were not reiined before being tested by Profiles since 
experience has shown that refinement does not substantially 
improve the profile of a loop Uiat is misfolded in the 
unrefined model. 

The disulphidc bonds in the cystine knot were formed 
between Cys26-Cys68, Cys57-Cysl02 and Cys60-Cysl04. A 
dimeric model was constructed in the same orientation as the 
PDGF dimer and intersubunitdisulphide bonds formed between 
Cys51 and Cys60. This seemed the most likely orientation 
because the Cys residues involved in the intersubunit disulphide 
bond are conserved in VEGF and PIGF and it has been shown 
by site-directed mutagenesis of these Cys residues that a similar 
pair of disulphide bonds is necessary for the dimerization and 
biological activity of VEGF (Potgens et al, 1994). 

The bond lengths and angles at the splice points where 
loops were added were refined using the Splice Repair facility 
of the Homology module. All residues except the two next to 
the splice points were fixed, a torsion force was placed on the 
0) angles of the peptide bonds at the splice point and 100 
steps of steepest descent minimization were performed. The 
conformations of mutated side chains were then adjusted using 
an automated procedure that minimizes the energy of side 
chain conformations by a systematic search of rotamer con- 
formations. 

A Profiles (Luthy et ai, 1992) analysis of the unrefined 
VEGF model showed that residues 17 and 18 were in very 
incompatible environments. These hydrophobic residues may 
interact with the hcparin-binding domain. These hydrophobic 
residues lie on tlie opposite face of the dimer to the receptor- 
binding loops and they may interact with the heparin-binding 
domain. Residues 17 and 18 were cut from the model and the 
dimer was re-refined. 

The dimeric model was solvated with a 9 A shell of T1P3P 
water (Jorgensen et at., 1983). Molecular mechanics and 
dynamics were performed using CHARMM with a united 
atom forcefield (Brooks et aL, 1983). The solvated model was 
energy minimized using the Adopted Basis Newton-Raphson 



algorithm with the CHARMM forcefield (Brooks et ah, 1983). 
Initially only hydrogens were minimized, then water molecules, 
followed by side chains and then the entire model. The model 
was then subjected to molecular dynamics. A 1 fs time step 
was used throughout the molecular dynamics simulation. The 
non-bonded cut-off distance was 8 A and the dielectric was 
distance dependent. The solvated model was heated from 0 to 
300 K over 1 ps and then equilibrated for 10 ps. This was 
followed by free dynamics at 300 K for 25 ps. Bonds involving 
hydrogen atoms were constrained. The model was then 
quenched from 300 to 0 K over 1 ps and reminimized without 
constraints. 

Stereochemical analysis of the models was performed using 
PROCHECIC (Laskowski et 1993). .The quality of the 
models was assessed using Profiles (Luthy et «/., 1992). 
Electrostatic potential grids were created using DelPhi (Sharp 
and Honig, 1990). All protein charges were taken from the 
AMBER charge templates provided with the DelPhi program 
(Sharp and Honig, 1990). 

Results and discussion 

MOLSCRIPT diagrams of VEGF and PIGF are shown in 
Figures I and 2. The dimer structure is similar to that in 
PDGF but is more curved along the long axis (Figures 3 and 
4). The r.m.s differences were calculated between the starting 
and refined structures (Table 1). These indicate that the refined 
structures have changed somewhat from the initial structure. 
The side chain nm.s. differences are large since the sequence 
identity between PDGF and VEGF/PIGF is low and the dimer 
lacks a hydrophobic core in which the conformational flexibility 
of the side chains would be restricted. A comparison of the 
structures of PDGF, TGFP2 and NGF shows that their back- 
bones diverge rapidly from each other away from the cystine 
knot. It seems reasonable, therefore, that the backbones of the 
refined VEGF/PIGF models diverge from that of the PDGF 
template. Ramachandran plots of the backbone angles (Figures 
5 and 6) in the models show that they fall into commonly 
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Hg. 2. MOLSCKIPT siereo drawing ol' die PlGF dimer. Hydrophobic residues in the dimer interiitce are drawn in ball-Hnd-stick representation. 
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l^1g. 3. Ribbon diograms of the V£GF dimer related by a 90" rotation, drawn using MOLSCRIPT. The cysteine residues that form the cystine knot and the 
intersubunit disulphide bonds are drawn us balls and sticks. 



observed regions of space. The conserved glycine in the 
cystine knot has a positive 0 conformation as is observed in 
PDGF and TGFP2. The disulphide bonds in the cysteine knot 
have similar conformations to those in PDGF. The p bulge 
preceding Cys97 in PDGF is also present before Cysl02 in 
VBGF and PIGR The dimer has the shape of a shallow trough 
with loop 1 projecting from the concave surface. The structure 
of loop I was built based on the structure of residues A39- 
A48 of PDB entry 3rp2. The systematic classification of (3 
hairpins (Sibanda et a}., 1989) predicts that the deletion of 
three residues in loop U will change the class 2:2i' hairpin in 
loop II of PDGF to a class 3:5 hairpin. Loop 11 was built using 
a class 3:5 hairpin from Isgt (Figure 7). The hydrogen bonding 
pattern in the P sheets is similar to that in PDGF. There are 
two hydrogen bonds between Lys32 and Met65 {Ser30 and 
Leu63 in PIGF), causing a bulge in the p strand following 
loop 1. Mutation of Met65 (Leu63 in PIGF) to a non- 



hydrophobic residue may change the positions of loops 1 and 
11 relative to each other since this residue is buried along 
the interface between the loops. The twofold symmetry of 
the dimer causes the clustering of the loops at the ends of the 
dimer. In the other cystine knot growth factors the residues 
involved in receptor binding are located in the loops. If diis 
is the case in PIGF and VEGF. then the loops will form two 
identical binding regions at opposite ends of the dimer that 
could bind two receptor molecules and induce receptor 
dimerization and activation. Mutations in loop 1 of PDGP do 
not disrupt the structure of the rest of the protein and this may 
also be the case in VEGF/PIGF, since this loop points into the 
solvent and does not interact with the body of the protein. 

It has been shown that VEGF is only active in the dimeric 
form and that dimerization is dependent on the formation of 
an inlersubunit disuIphide.Jbond between Cys51 of one subunit 
and Cys60 of the other (Potgens et al, 1994). Mutants lacking 
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Fig. 4. Ribbon ditignifns uF the PICF dtmer related by h 90' rotation, drawn using MOLSCRIPT. Thu cysteine residues that fonn the cystine knot nnd the 
iniereubunil disulphide bunds iire dniwn as bulls and sticks. 



Ramachandran Plot 
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Fig. 5. Kumachandron plot for VEGF. prepared using PROCHBCK. 

either of these' Cys are unable to fonn dimers. The coexprcssion 
of CysSl-^SerSl mutants with Cys60->Ser60 mutants leads 
to formation of heterodimers that are fully active despite 
having only one intersubimit disutphide bond. 

The cysteine knot is formed by the disulphide bonds Cys26- 
Cys68, Cys57-Cysl02 and Cys61-Cysl04. These disulphide 
bonds have similar conformations to those in PDGF. The 
mutation of any one of Cys26, Cys57, Cys6l or Cys68 to a 
Ser residue results in impaired dimerization (Potgens et al, 
1994; Claffey et al, 1995). This cannot be due to disruption 
of an intersubunit disulphide bond involving these Cys residues 
since they are not in a position to form disulphide bonds with 
any of the Cys residues in the symmetry-related subunit. These 
mutations disrupt the structure of the cystine knot, causing a 
change in the position of Cys60. which is located within the 
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Fig. 6. Romscliandmn plot for PlGF, prepared using PROCHECK. 

backbone ring of the knot. Cys26 forms one half of the 
disulphide bond that penetrates the ring of the cystine knot; 
mutation of Cys26 may disrupt the tight packing in the knot. 
One study (Potgens et ah, 1994) found that the mutant 
Cys57->Ser57 can still form dimers at the level of 50% of the 
wild type and that these dimers are biologically active but 
another (Claffey et ai, 1995) found that, like the other mutants, 
this mutant is imable even to dimerize. Cys57 is at the end of 
P strand 2 next to loop 11 and may be less crucial to the 
integrity of the fold; mutants at this site might still be correctly 
folded and able to form noncovalent dimers. It is interesting 
that Claffey et al (1995) found that the mutation Cysl02->Ser 
prevents the synthesis and/or secretion of this isoform. In 
addition, this isoform is a dominant-negative mutant that limits 
the synthesis and/or secretion of wild-type VEGF. This suggests 
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Klg. 7. Buckbono stnicluro urihc loop frum PDB cnlry Isgt thut was used 
lu model loop 111. Douud lines indicalu hydrogen bonds. 



(A) 



VEGF monomer 




Sequence number 



(B) 



VEGF dimer 




Sequence number 

Fig. 8. Profile plots for the monomeric (A) and dimeric (B) fonns of VEGF. 
A seU-compatibility score {S) is determined Tor each residue in the 
sequence. The compatibility of the sequence with the structure can be 
assessed by plotting the .V score against residue number In a window. 
Negative iV scores indicate putQtively misfolded regions. A smoothing 
window size of 21 residues is used end .scores for the Hrst ten and last ten 
residues have no significance. 

that the Cysl02->Ser mutant can form a heterodimer with 
normal VEGF that is subsequently degraded intracellularly. 
The stability of the dimer may also be affected by inter- 
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Fig. 9. Profile plots for the monomeric (A) and dimeric (B) forms of PIGF. 
Sec the legend to Figure 8 for an explanation of the ploL 

subunit hydrogen bonding. There are several charge pairs 
buried along the dimer interface: in VEGF Lys48 forms charge 
pairs with Asp63 of the other monomer while Arg82 forms a 
charge pair with Glu64. in PIGF there are charge pairs formed 
by Lys82 with Glu64 and GluBO with Arg27. A striking feature 
of the models is that almost none of the hydrophobic residues 
have more than 90% of their accessible surface area buried in 
the interior of the protein. The only exception is Val52. which is 
completely buried (0% of surface area accessible to solvent) 
in both VEGF and PIGR This is adjacent to Cys51 and 
the burial of this residue may play a part in determining the 
conformation of the cystine knot. However, dimer formation 
does bury some hydrophobic surface area, as can be seen from 
a comparison of the profiles of the monomeric and dimeric 
forms of VEGF and PIGF (Figures 8 and 9). Hydrophobic 
residues involved in the dimer interface are listed in Table II. 
The residues at positions 20, 29, 43, 46, 52, 78, 81, 83 and 
96 are involved in the dimer interface in both proteins. A 
model of the VEGF/PIGF heterodimer might provide further 
insights into the requirements for dimerization. The models 
show that the PDGF-like domains of VEGF and PIGF are 
essentially two-dimensional structures that lack a hydro- 
phobic core. 

The Profiles suite of programs (Luthy et al., 1992) was used 
to assess the quality of the models. Profiles measures the 
compatibility of each residue in a model with its environment 
and calculates a Z score that is a measure of the compatibility 
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fig. 10. Solvenl-accessible surfhces of the VEGF (A) and PIGP (B) monomers coloured according; to the electrostatic energy at each point on the surface. 
The energies wen; taken from an cicclrustatic grid atlculated using the DulFhi module in tnsightU. The units of the spectrum are k'l\ 



of the entity sequence with the fold of the model. Correctly 
folded proteins have a Z score at least 45% of the maximum 
theoretically possible score ((he S score). The dimeric VBGF 
model has o 2 score that is 57% of the S score. The residue 
by residue profile (Figures 8 and 9) does not indicate any 
misfolded regions although the scores for loops i and ill are 
low. The dimeric PIGF model has a Z score that is 70% of 
the ^ score. While it is possible that the N-terminal region of 
the domain that could not be modelled interacts with the part 
modelled here, the Z scores and the profiles indicate that the 
fold is basically correct. The areas of greatest error are likely 
to be the loops, especially the large loop 1. A comparison of 
the profiles of the dimeric and monomeric models (Figures 8 
and 9) shows that dimer formation improves the scores 
along the dimeric interface by burying hydrophobic residues. 
Dimerization buries about 900 of surface area per monomer. 
Of this, about 380 A^is buried along the dimer interface; this 
is less than in PDGF. Unlike PDGF. VEGF is unable to forni 
stable, biologically active noncovalent dimers. This may be 
because the smaller surface area buried in VEGF is insufficient 
to .stabilize the dimer in the absence of intersubunit disulphide 



bonds. By contrast, the NGF dimer, which Is non-covalent, 
buries 1166 A- of surface area per monomer. 

The receptor-binding domains of VEGF206 have been studied 
using monoclonal antibodies raised against peptides from the 
sequence of VEGF206 (Sioussat e( aL, 1993). Using the Chou 
and Fasman (1978) algorithm, peptides were predicted to 
contain elements of regular secondary structure. The study 
identified the N- and C-lermini of VEGF206 as being important 
in receptor binding. If the models presented are correct then 
none of the peptides used would include complete loops and 
the |5 strands joined by the loops. Since antibodies to proteins 
are generally formed against three-dimensional structure rather 
than sequence, antibodies raised against the peptides used in 
diis study would not recognize the loop regions. None of the 
available secondary structure predictors is able to predict 
accurately the secondary structure of the cystine knot growth 
factors and are therefore not useful guides to designing peptides 
for a study of this kind. 

Delphi calculations of the electrostatic potential over the 
surface of the dimers reveal a number of differences in the 
distribution of positive and negative charges (Figure ID). Loop 
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His 86 



VEGF monomer 

fig. 11. MOLSCRIPT ribbon diagram of the VEGF monomer Residues in 
the loops that are postulated to be involved in receptor binding are drawn In 
ball-and-stick represeniaiion. 




PIGF monomer 

Hg. ]2. MOLSCRIPT ribbon drawing of the J*1GF monomer with residueit 
postuluted lo be involved in receptor binding drawn In ball-and-stick 
representation. 

1 in VEGF contains a large, elongated negatively charged 
region formed by Asp34, Gln37, Glu38, Pro40, Asp41, Glu42 
and Glu44, whereas this loop in PIGF is almost uncharged. 
The body of the PIGF dimer contains a large area of positive 
charge consisting of Arg55, His99, Arg31, His67, Pro70 and 
Ser97, which is absent in VEGF The residues Lys84, His86 
and l-iis90 in VEGF form an elongated positively charged 
patch at the tip of loop 111. The corresponding hairpin in PIGF 
has a positively charged surface formed by the side chains of 
Arg84 and Arg88 (Figures 10 and 12); this surface is smaller 



than that in VEGF. This loop also includes a negatively charged 
side chain, Asp87; the corresponding residue in VEGF is Gin. 
In both VEGF and PIGF, loop 11 has a negatively charged 
patch formed by Asp63. Glu64 and Glu67 in VEGF and by 
Asp63 and Glu64 in PIGR 

There are several patches of exposed hydrophobic residues 
on tlie surfaces of VEGF and PIGF. The residues MetSl. i)e83 
and !le91 form a hydrophobic patch on the surface of loop HI 
in VEGF. tn PIGF the residue corresponding to Met8I (Leu8I) 
is buried and Ile91 is replaced by a Tyr residue, resulting in 
a much smaller hydrophobic patch formed by lle83. The 
conserved Leu66 is exposed in both VEGF and PIGF. There 
is a large area of hydrophobic surface formed at the interface 
of loops I and II by Val33, Val35, Phe37, Leu80 and Val92 in 
PIGF and by Val33, Ile35. Phe37 and lleSO in VEGF 

The charged residues in VEGF that are involved in receptor 
binding have been identified by alanine-scanning mutagenesis. 
Arg82, Lys84 and His86 in VEGF are critical for binding to 
KDR/fllC'I while Asp63, Glu64 and Glu67 were implicated in 
binding to flt-1 and to a lesser extent, KDR (Keyt et ai, 1996). 
A model of VEGF was constructed by aligning the sequences 
of VEGF and PDGF and substituting the side chains of these 
charged residues for the corresponding residues in PDGF. This 
model does not contain coordinates for loop 1 and was not 
subjected to refinement by ene(;gy minimization or molecular 
dynamics. The backbone structure of this model is identical 
with that of PDGF. The model of VEGF presented here shows 
that the backbones of VEGF and PDGF diverge significantly 
from each other outside the regioii of the cystine knot. 

This model locates the positively chaiiged residues Arg82, 
Lys84 and His86 in a cluster on the edge of the loop 111 |3 
hairpin. The model presented here shows that the side chains 
of Arg82, Lys84, His86 and His90 form a positively charged 
surface on the plane of the P hairpin, with His86 forming part 
of the P hairpin (Figures 10 and II). Although His90 forms 
part of this positively charged surface mutation of this residue 
to Ala has little effect on binding. The position of the negatively 
char^ged residues Asp63, Glu64 and Glu67 is similar in both 
models. The potential extra glycosylation site introduced by 
mutation of Glu64 is located at the tip of loop II, which mediates 
binding to flt-L Glu64 is one of the residues implicated in 
binding to fit-l\ it is therefore not surprising that mutation of 
this residue to Asn reduces binding to flt-1 -40-fold. No effect 
on binding lo either^?-/ or ICDR was detected by mutagenesis 
of charged residues in loop 1. The comparison of the models 
of VEGF and PIGF suggests that the inability of PIGF to bind 
K.DR may be due to differences in the positively charged 
surface in loop HI of PIGF as compared with VEGF and the 
absence of exposed hydrophobic residues in loops I and II. 
The ability of PIGF to bind flt-l is consistent with the presence 
of a negatively chaiiged surface in loop 11 similar to the one 
that mediates VEGF binding to flt-L 

VEGF has a putative N-linked glycosylation site with the 
consensus sequence Asn-Thr-Leu, at residues 15-11. Mutation 
of Asn75 to Tyr alters the electrophoretic mobility of the 
VEGF dimer, consistent with a lack of glycosylation (Claffey 
et al, 1995). The VEGF model shows that the side chain of 
Asn75 is exposed to solvent and is located near the middle of 
the fourth p strand (Figure II). Mutation of this site has no 
effect on receptor-binding activity (Claffey et aL, 1995), which 
is not surprising given that it is not located near the receptor- 
binding loops. The site is conserved in PIGF. 

It is hoped that the models presented here will be useful 
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Computer models of VEGK and PIGF 
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til II 
SIEEAVP AVCKTRTVIYEIPR SQVDPTSANF LIWPPCVEVKRCTGCC 
lAEPjijyil AECKTRTEVFEISR RLIDRTNANF LVWPPCVEVQRCSGCC 



FMDVYQRSEffiEimimiZQEYPD- 
FMDVYQRS YCRP lETLVD IFQ EYPD- 
FMDVYORS YCHP lETLVD IF OEYPD- 
FQEVWGRSYC MX^ERXiVDVV SEYPS- 
FNEVWGRS YCRPMEKLVYIA DEHPN- 



- E I EF I FKP S CVP LMRCG GCC 
■EIE YIFKPSCVPLMRCGGCC 
■ E I E Y I FKP S CVg l^mQhGQQ 
■EVEHMFSPSCVSLLRCTgCG 
-EVSHIFSPSCVLLSRCSGCC 
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NTSS VKCOPSRVHHRSVKVAKVEY VRK KPKLKEVOVRLEEHLECACATTSLN 

NNRNVQCRPTQVQJiRFVQVRKIEJiVRKKPIFKKATVTI^EDHLACKCBTVM- 

NDEG LECVPTEESNITMOIMRIK PHO G QHIGEMSFLQHNKCECRPKKD R 

NDESL ECVPTEEFNITMOIMRIK PHO— S OHIGEMSFLOHNKCECRPKKD K 
NDEA LECVPTSESNITMOIMRIK PHO— S OHIGEMSFLQHSRCECRPKKDR 
GDENLHCVPVfiTAMVTMPLLKIRSGD"--RPSYVELTFSQlHVRCECRgLREK 
GDEGLHCVALKTANITMOILKIPPNR— DPHSYVEMTFSQDVLCECRPILET 
. * it 1r 



Fig. 13. Clustal sequence alignment ol'the sequences oF human PDGF-A and PDGF-B, human, bovine and mouse VEGF and human and rot forms of PIGF. 
The alignment begins at residue 19 of the mature VEGF sequence and ends at residue 109 of that sequence. This is the region of VEGF corresponding to the 
model. Cysteine residues involved in the cystine Icnot arc in bold type. Conserved residues are marked by an asterisk and conservative substitutions by o 
period (.)• Residues underiined ore in the p sheets. 



Tabic I. Root mean square ditTerences between storting and refined 
structures of PIGF and VEGF 





Atom selection 


R.m.s.d. (A) 


PIGF 


All non-hydrogen atoms 


3.53 




Backbone atoms 


1.64 




oc-Carbon atoms 


1.60 


VEGF 


All non-hydrogen atoms 


2.13 




Buckbunu Htoms 


1.90 




a-Carbun atoms 


1.82 



Table II. Hydrophobic residues with surface urea buried in the dimer 
interfuses of VEGF and PIGF 



VEGF 


PIGF 


vui2a 


Val20 


ile29 


Leu29 


Leu32 


Uu32 


Phe36 




Ile43 


Val43 


VatS2 


Vol52 


Leu54 


LeuS4 




LeuSS 


Met78 


Met78 


lleKO 




Met81 


Uu81 


Ile83 


Ile83 


Iie91 




Mcl94 




Phc96 


Phc96 



as guides to the design of future mutagenesis cxpeiiments. 
Modelling of the receptors is in progress. In combination with 
models of the receptors, the models of VEGF and PIGF will 



be used to investigate the determinants of receptor binding. 
An understanding of these would be a first step towards the 
rational design of VEGF inhibitors which might prove useful 
as inhibitors of angiogenesis and tumour growth 
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